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"Strain engineering" in functional materials has been widely explored to tailor the physical properties of 
electronic materials and improve their electrical and/or optical properties. Here, we exploit both in plane 
and out of plane uniaxial tensile strains in MoS 2 to modulate its band gap and engineer its optical properties. 
We utilize X-ray diffraction and cross -sectional transmission electron microscopy to quantify the strains in 
the as-synthesized MoS 2 nanosheets and apply measured shifts of Raman-active modes to confirm lattice 
strain modification of both the out-of-plane and in-plane phonon vibrations of the MoS 2 nanosheets. The 
induced band gap evolution due to in-plane and out-of-plane tensile stresses is validated by 
photoluminescence (PL) measurements, promising a potential route for unprecedented manipulation of the 
physical, electrical and optical properties of MoS 2 . 



Transition metal dichalcogenides, such as MoS 2 , have attracted considerable attention for their potential 
advantages in catalysis 14 , transistors 5 7 , batteries 8,5 , gas sensors 10 , photodetectors 11 , and other optoelectronic 
devices 12 . Bulk crystals of MoS 2 , with an indirect band-gap of 1.3 eV, possess strong intra-layer sulfur-metal 
covalent bonds where the metal atoms are sandwiched by sulfur atoms but their interlayer stacking is weakly 
supported by van der Waals forces 13 . Single-layer MoS 2 , however, shows distinctly promising semiconductive 
properties with a direct band-gap of 1.9 eV 14 , exhibiting strong photoluminescence 15 , high in-plane carrier 
mobility (200-500 cm 2 /V - s) 16 and robust mechanical properties 17 . Interestingly, such indirect-to-direct band 
gap transition from the bulk form to the monolayer form due to quantum confinement leads to strong enhance- 
ment of photoluminescence (PL) from MoS 2 14 , therefore creating a series of promising applications for two- 
dimensional (2D) material-based optoelectronics 1118 . This is afforded through the large yield of monolayer MoS 2 
synthesis by chemical vapor deposition (CVD) that has been utilized to build up atomically thin optoelectronic 
devices 19,20 . Moreover, it would be desirable to achieve a wider optical spectral response if the band gap of layered 
MoS 2 can be systematically controlled. "Strain engineering" has been shown to exert a profound impact on the 
properties of MoS 2 21 ~ 23 . Since MoS 2 can endure large strains before breaking 24,25 , tuning the MoS 2 band gap by 
lattice strain becomes an important strategy to enhance the performance of nanodevices made from MoS 2 . In this 
work, we systematically investigated the effects of in-plane and out-of-plane tensile strains of MoS 2 nanosheets 
and explored the associated strain engineering of their band gap. The lattice strains in as-synthesized MoS 2 
nanosheets were quantified using the X-ray diffraction (XRD), cross-sectional high-resolution transmission 
electron microscope (TEM) analysis and Raman spectroscopy. The band gap evolution of MoS 2 nanosheets 
induced by the strains was validated by PL measurements, which is qualitatively consistent with our theoretical 
analysis using density functional theory (DFT) calculations. 

The phase of MoS 2 nanosheets synthesized using a hydrothermal method (Fig. SI) is confirmed to be a 
hexagonal structure of MoS 2 with space group P 63/mmc as indexed by JCPDF card number 37-1492. 
Particularly, the peak of (002) plane (c-axis) shifts to a smaller angle while the peak of (100) plane shifts to a 
larger angle, indicating that there is an expansion along [001] direction and in-plane compression in the as- 
synthesized MoS 2 nanosheets. The atomic plane spacing in the MoS 2 nanosheets can be calculated from the XRD 
pattern using the Bragg equation 26 : 

2dsin9 = nX (l) 
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where d is the atomic plane spacing, 0 is the Bragg scattering angle, n 
is the diffraction order and A is the wavelength of incident X-ray. 
Using Equation (1), the XRD result showed a calculated interlayer 
spacing ((002) Miller indices spacing) of 6.1924 A for the as-synthe- 
sized MoS 2 nanosheets that is larger compared to that of bulk MoS 2 
(6.1554 A) 27 . The observed difference in the (002) plane spacing 
causes a uniaxial tensile strain along [001] direction in the as-syn- 
thesized MoS 2 nanosheets. By tuning the PH value in the reaction 
solution, the lattice strains in as-synthesized MoS 2 nanosheets could 
be modified (Supporting Information, Figure S2). The morphology 
of the as -synthesized MoS 2 nanosheets was characterized by scan- 
ning electron microscopy (SEM), showing that the 2D MoS 2 
nanosheets rolled up and thus formed a three-dimensional (3D) 
network with highly curved walls (Figure S3). From the high resolu- 
tion TEM (HRTEM) image of the basal plane of the as-synthesized 
MoS 2 nanosheets shown in Figure lb, we can see hexagonal rings of 
alternative molybdenum and sulphur atoms in each unit, meaning 
that the MoS 2 nanosheets are of high quality with defect- free atomic 
lattices. The corresponding zone axis is along [001] direction as 
confirmed using the fast Fourier transform (FFT) shown in the inset 
of Figure lb. The marked (100) and (010) atomic plane distances 
(~2.7 A) are both smaller (—1%) than their corresponding lattice 
spacing in bulk MoS 2 28 , indicating an in-plane biaxial compressive 
strain in the basal plane. We also applied HRTEM to image the cross- 
section of the MoS 2 nanosheets (Figure lc). Typical nanosheets con- 



sist of 5, 6, 9, and 10 layers, from which the interlayer spacing can be 
measured accordingly (Supporting Information). The average of the 
measured interlayer spacing is 6.2282 A in the as-synthesized MoS 2 
nanosheets (Figure Id), which is about 1% larger than the bulk inter- 
layer spacing (6.1554 A) indicating an out-of-plane uniaxial tensile 
strain along [001] direction, consistent with the XRD results dis- 
cussed above. These results confirm that the as-synthesized MoS 2 
nanosheets contain an out-of-plane uniaxial tensile strain and an in- 
plane biaxial compressive strain as depicted in Figure le. 

Figure 2a shows Raman spectra of as-synthesized MoS 2 nano- 
sheets (grown from pH = 7 and pH < 7 solutions, respectively), 
strain-partially released MoS 2 nanosheets (grown at pH < 7 and 
transferred by solution method to Si0 2 /Si substrate, see supporting 
information), bulk MoS 2 (SPI Supplies) and monolayer MoS 2 (CVD 
growth on Si0 2 /Si substrate). Compared with bulk MoS 2 , there is a 
red shift of A lg mode observed in the monolayer MoS 2 with Aro = 
3 cm -1 , and for E 2g mode, a blue shift (Aco = 2 cm -1 ) exists in the 
monolayer MoS 2 . In the Raman spectra of MoS 2 nanosheets grown 
from pH = 7 solution and PH < 7 solution, only red shifts for both 
Ejg and A lg modes were observed, compared to bulk MoS 2 . Com- 
pared to as-synthesized MoS 2 nanosheets, blue shifts of both E 2g and 
A lg modes in the strain-released MoS 2 nanosheets were observed. 
Figure 2b shows the frequencies of E^ and A lg modes of MoS 2 
nanosheets extracted from Figure 2a. The largest red shifts of both 
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Figure 1 (a) XRD pattern of MoS 2 nanosheets grown from PH < 7 solution with the line indexed by JCPDF card number 37- 1492. (b) HRTEM image of 
the basal plane of MoS 2 nanosheets grown from PH < 7 solution. Atomic lattice is illustrated by the hexagonal rings of alternating molybdenum and 
sulphur atomic sites in each unit as illustrated by the cartoon color spheres. The inset shows the corresponding fast Fourier transform (FFT) with a zone 
axis of [001]. (c) Cross-section HRTEM images of the MoS 2 nanosheets grown from PH < 7 solution. Each layer is denoted by the yellow spheres 
representing the Mo atoms. The direction of the cross-section HRTEM is along [001]. (d) The average of measured interlayer spacing from HRTEM 
pictures for different MoS 2 nanosheets consisting of (a) 5, (b) 6, (c) 9 and (d) 10 layers, respectively. The solid line represents the bulk interlayer spacing, 
(e) Diagram of in-plane biaxial compressive strain and out-of-plane tensile strain in the as-synthesized MoS 2 nanosheets grown from PH < 7 solution. 
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Figure 2 (a) Raman spectra of as-synthesized MoS 2 nanosheets including grown from pH = 7 (in pink) and pH < 7 (in blue) solution, strain-partially- 
released MoS 2 nanosheets (in gray), bulk MoS 2 (SPI Supplies) (in red) and monolayer MoS 2 (CVD growth on Si0 2 /Si substrate) (in black). 
The inset images near each spectra correspond to their sample morphologies. Scale bar is 300 nm. (b) The frequencies of E 2g and A lg modes extracted 
from (a). The Raman peak shift error bar indicates the spectrometer resolution, (c) PL spectra of as-synthesized MoS 2 nanosheets (grown from 
pH = 7 and pH < 7 solution), strain-partially-released MoS 2 nanosheets, bulk MoS 2 (SPI Supplies) and monolayer MoS 2 (CVD growth on Si0 2 /Si 
substrate), (d) The band-gap shifts induced by different lattice strains in MoS 2 . 



E 2g and A lg modes are observed in the MoS 2 nanosheets grown from 
PH = 7 solution. The strain-caused band-gap shift was also validated 
as shown in Figure 2c where PL spectra are plotted for bulk MoS 2 , 
monolayer MoS 2 , strain-partially released MoS 2 nanosheets and as- 
synthesized MoS 2 nanosheets. From the PL spectra of the as-synthe- 
sized MoS 2 nanosheets, two peaks at 693.9 nm (1.787 eV) (in blue) 
and 684 nm (1.813 eV) (in pink) were clearly observed, correspond- 
ing to a direct transition between the top valence band K point and 
the bottom conduction K points in the Brillouin zone 14 . In compar- 
ison with MoS 2 nanosheets grown from PH < 7 solution, a blue shift 
(26 meV) of the direct band transition was observed in MoS 2 
nanosheets grown from PH = 7 solution synthesized MoS 2 
nanosheets, while a red shift (21 meV) of the direct band transition 
was observed in the strain-partially released MoS 2 nanosheets. In 
monolayer MoS 2 , most of the excited electrons and holes recombine 
with each other through direct band-to-band transition at the K 
point in the Brillouin zone, leading to strong photoluminescence 
near 1.87 eV. A weak feature was also observed at 613 nm 
(2.023 eV as noted) due to direct-gap hot luminescence 14 . For bulk 
MoS 2 , an indirect band gap transition at 924.7 nm (1.341 eV) was 
observed, which is a transition between the top valence band T point 
and the bottom conduction band halfway between T and K points in 
the Brillouin zone 29 . The band-gap shifts due to different strain 



effects in MoS 2 that were extracted from Figure 2c are shown in 
Figure 2d. 

To further illustrate the band-gap modulation by strain, we inves- 
tigated the influence of in-plane uniaxial tensile strain on the phonon 
spectra and the band gap of exfoliated MoS 2 nanosheets. The exfo- 
liated MoS 2 nanosheets were transferred onto a flexible substrate 
(details in Supporting Information). The solidification process for 
the polydimethylsiloxane (PDMS) flexible substrate caused the exfo- 
liated MoS 2 nanosheets to strongly adhere to the surface of PDMS 
(Figure S7). Upon bending the PDMS substrate, an in-plane uniaxial 
tensile strain is exerted at the top-most plane containing the exfo- 
liated MoS 2 nanosheets. As shown in Figure 3a, the applied strain 
significantly causes blue shifts in A lg and E 2g modes of the exfoliated 
MoS 2 nanosheets. Assuming that the bending induced strain is 
entirely transferred to exfoliated MoS 2 nanosheets, we estimated 
the in plane uniaxial tensile strain s existing in the exfoliated MoS 2 
nanosheets by calculating the ratio of the PDMS substrate thickness d 
to twice the radius of the curvature R as 8 = d/2R 30 . Figure 3b shows 
the frequencies of E 2g and A lg modes of the exfoliated MoS 2 
nanosheets with different local strains extracted from Figure 3a. In 
Figure 3c prior to bending the PDMS substrate, the exfoliated MoS 2 
nanosheets exhibit a PL peak at 880.1 nm (1.409 eV), which is sim- 
ilar to the indirect band transition in bulk MoS 2 . We also observed 
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Figure 3 | (a) Raman spectra of the exfoliated MoS 2 nanosheets on PDMS flexible substrate. The insets in (a) show the schematic illustrations of the 
exfoliated MoS 2 nanosheets on flexible substrate undergoing strain effect. The in-plane uniaxial tensile strain was introduced by bending the PDMS 
flexible substrate, (b) The frequencies of E 2g and A lg modes extracted from (a). The Raman peak shift error bar indicates the spectrometer resolution, 
(c) The PL spectra of exfoliated MoS 2 nanosheets with 0% (in black), 0.47% (in red), 1.21% (in blue) in-plane uniaxial tensile strain, (d) The PL peaks of 
exfoliated MoS 2 nanosheets on flexible substrate extracted from (c). The band-gap is decreases by an increase of in-plane uniaxial tensile strain. 



another direct band transition PL peak at 695.1 nm (1.784 eV), cor- 
responding to the same direct band-to-band transition as observed in 
as-synthesized MoS 2 nanosheets. With the in-plane uniaxial tensile 
strain, red shifts of direct band-to-band transition and indirect band 
transition are observed by bending the PMDS substrate. Figure 3d 
shows the PL peaks of the exfoliated MoS 2 nanosheets with different 
local strains extracted from Figure 3c. 

MoS 2 consists of weak bonds between different S-Mo-S layers as 
illustrated in the stacked structure of Figure 4a. Each layer is elec- 
trically neutral, such that nature of the forces between adjacent layers 
is van der Waals' type 31 . The van der Waals forces between S atoms in 
adjacent layers dominates the interlayer forces. Figure 4b shows 
atomic displacements of the four Raman-active modes. In the out- 
of-plane A lg mode, the restoring force is primarily due to interlayer 
van der Waals interaction 32 . With increasing the number of layers, 
the additional "spring" between S atoms in neighboring layers 
enhances the restoring force, resulting in an increase of the A lg mode 
frequency. The in-plane E 2g mode is a symmetric mode, which 
vibrates out of phase as sketched in Figure 4b. The presence of the 
interlayer interactions can also increase the effective restoring force 
acting on the basal plane atoms 31 " 33 . It has been demonstrated that in 
bulk graphite, this kind of inter-layer interactions affect the intra- 
layer bonding and stiffen the in-plane lattice vibrations 34 . Similarly, 
with an increase of number of layers in MoS 2 there is an increase in 



the frequency of the in-plane E 2g mode. Compared to the Raman 
spectrum of bulk MoS 2 , red shifts are observed for both A lg and E 2g 
modes of the as-synthesized MoS 2 nanosheets. As characterized 
above, there is a uniaxial tensile strain along [001] direction in the 
as-synthesized MoS 2 nanosheets grown from PH < 7 solution. The 
out of plane uniaxial tensile strain weakens the interlayer interac- 
tions, leading to a decrease in frequencies of both A lg and E 2g modes. 
The inset images in Figure 2a show that the average pore size of 3D 
MoS 2 networks grown from PH = 7 solution is —100 nm, smaller 
than in 3D MoS 2 networks grown from PH < 7 solution. In terms of 
lattice strain, the difference in the curvature causes a larger strain 
effect in MoS 2 nanosheets grown from PH = 7 solution as supported 
by XRD (details in Supporting Information) and the red shifts of 
both A lg and E 2g modes observed in the Raman spectra, compared to 
that of MoS 2 nanosheets grown from PH < 7 solution. During the 
sonication process, the 3D networks can be cracked to smaller sheets 
of strain-released MoS 2 as indicated in the inset of Figure 2a, causing 
blue shifts of A lg and E 2g modes compared to as-synthesized MoS 2 
nanosheets. However, mild red shifts (1 cm -1 ) of both A lg and E 2g 
modes were still observed in strain-partially released MoS 2 
nanosheets, in contrast to bulk MoS 2 , due to residual strains in these 
discontinuous nanosheets. This residual strain resulted from the 
small curvatures observed in the sonicated MoS 2 nanosheets with a 
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Figure 4 | (a) Crystal structure of hexagonal MoS 2 ,with strong intra-layer sulfur-metal covalent bonds and weak inter-layer van der Waals stacking, 
(b) Illustration of the atomic displacements of the four Raman-active modes, (c) The calculated band-gap energies for MoS 2 nanosheets with a different 
number of layers modulated with 1% out-of-plane tensile strain, (d) The simulated Raman scattering of bulk. MoS 2 with the effect of 1% out of plane 
tensile strain and 1 % in-plane tensile strain. 



morphology similar to broken-bamboos (Figure S8, Supporting 
Information). Note that the increase in frequency of the A lg mode 
from the monolayer to the bulk MoS 2 is 3 cm -1 . One might expect 
that the additional interlayer van der Waals interactions should raise 
the frequency of E 2g mode with an increase of number of layers. 
However, according to our experimental results shown in 
Figure 2a, the E 2g mode in bulk MoS 2 has a lower frequency 
(~2 cm"') than in monolayer MoS 2 . This frequency decrease of 
the E 2g mode is possibly due to the long range Coulomb interac- 
tions 3235 from the coupled dipoles induced from Mo-S bonds. The 
Coulomb interactions tend to decrease the frequencies of both E 2g 
and E lu (TO) modes 32,35 . An increase in the number of layers leads to 
a stronger long-range Coulomb forces and stronger interlayer inter- 
actions that are competitive with each other. In few-layer MoS 2 
samples, these Coulomb forces are larger than the total forces arising 
from the interlayer interactions 35 , while in thicker layered samples, 
the interlayer interactions become stronger and surpass the Coulomb 
interactions. In monolayer or few layered samples, the Coulomb 
interactions dominate the Raman peak shifts. In this case, an in- 
plane tensile strain effect can lead to a red shift in Raman peaks as 
reported by Conley et al 21 , resulting from the increase of the Coulomb 
interactions. While in the thicker layered samples, the interlayer 
interactions become dominated for the shift of Raman-active mode 
frequencies. When the number of layers reaches up to 6, the layered 
sample retains bulk-like behavior and there is no continuous fre- 
quency decrease in the observed E 2g mode 14,34 . In our synthesized 
MoS 2 nanosheets, the Coulomb interactions are smaller than the 
interlayer interactions. Therefore, undergoing out of plane uniaxial 
tensile strain, the shifts in frequencies of A lg and E 2g modes are 



consistent with the classical model of coupled harmonic oscillators 33 . 
In addition, the larger in-plane tensile strain gives rise to larger blue 
shifts of both E 2g and A lg modes as demonstrated in Figure 3a and 3b. 
The in-plane uniaxial tensile strain can enhance the interlayer inter- 
actions in the exfoliated MoS 2 nanosheets, leading to stronger in- 
plane and out-of-plane effective restoring forces acting on the 
atoms 31,32 . Therefore, higher frequencies of A lg and E 2g modes are 
observed in the exfoliated MoS 2 nanosheets with an in-plane uniaxial 
tensile strain. 

It has been established that the energy states of the valence band 
and conduction band mainly originate from the 3p orbital of S atoms 
and 4d orbital of Mo atoms 36 . The interatomic distance change 
between the Mo and S atoms due to the strain effect can lead to a 
different superposition of those atomic orbitals. The variations in the 
overlap between the 3p z orbital of S and 4d z 2 orbital of Mo cause 
shifts in the energy states of the valence band and conduction band 36 . 
The presence of the lattice strains therefore modifies the overlap 
between the S orbital and the Mo orbital, which results in an energy 
shift in the band gap. In the presence of in-plane uniaxial tensile 
strain, the direct band gap is tailored by a rate of —35 meV per 
1%, while the indirect band gap is less sensitive to this tensile strain 
(Figure 3d). The strain is employed to tailor the band-gap by modi- 
fying the curvature of the bands, reflecting in a change in the effective 
masses of electrons and holes 37 . For MoS 2 , the location shift of con- 
duction band extremum generally occurs for much smaller strain 
than the location shift of the valence band extremum 37 . Therefore, 
the strain-induced change in the effective mass of the electron in the 
conduction band extremum halfway between F and K points in the 
Brillouin zone can reflect the indirect bang-gap evolution of MoS 2 . 
As subjected to in-plane uniaxial tensile strain, the effective mass of 
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electrons in the conduction band extremum halfway between T and 
K points exhibits a slow change with this tensile strain 37 . This eluci- 
dates that the MoS 2 indirect band-gap is less sensitive to in-plane 
uniaxial tensile strain. However, the effective mass of electrons at the 
bottom conduction band at K point decreases monotonically by 
increasing in-plane tensile strain. Therefore, the tunability of the 
indirect band gap is much smaller than that of the direct band gap 
in MoS 2 nanosheets by in-plane uniaxial tensile strain. 

The evolution of the band gap of the as-synthesized MoS 2 with the 
number of layers varied as 5, 6, 9, and 10 is also numerically analyzed 
with the influence of 1% lattice strain as shown in Figure 4c. The 
calculated band gap tunability by strain is —85 meV per 1% strain. 
The deviation between the experimental and numerical results is 
because the strain in the MoS 2 nanosheets is only partially released 
after the sonication process. Note that the calculated direct band gap 
value is smaller compared to the results in the PL spectrum. This is 
mainly due to that the standard exchange-correlation functions of 
the generalized gradient approximation (GGA) are used in the cal- 
culations, by which we usually obtain underestimated results for 
band-gap energies 33,38 . However, the band-gap computation indi- 
cates that the existence of the out-of-plane uniaxial tensile strain 
enhances the band gap of the MoS 2 , consistent with the experimental 
results. The simulated Raman scattering (Figure 4d) shows red shifts 
of Raman-active modes (A lg and E^) with a 1% out of plane tensile 
strain and blue shifts of Raman- active modes (A lg and E\) with a 1% 
in-plane tensile strain, which is also consistent with our experimental 
results. The decrease in PL intensity and peak broadening in as- 
synthesized MoS 2 nanosheets (Figure 2c) could be due to the expo- 
nential enhancement in non-radiative electron — hole recombination 
processes, reducing the probability of radiative transition 39 . 

In summary, we have experimentally demonstrated the lattice 
strain effects on the evolution of the band gap in as synthesized 
and mechanically exfoliated MoS 2 nanosheets. Concurrent studies 
of X-ray diffraction, cross-sectional high-resolution transmission 
electron microscopy analysis and Raman spectroscopy reveal the 
correlation between the structures of the MoS 2 nanosheets and their 
optical properties. The out-of-plane uniaxial tensile strain gives rise 
to a blue shift of MoS 2 band-gap whereas the in-plane uniaxial tensile 
strain results in a red shift of its band-gap. Our results suggest that 
manipulating the lattice strain by synthetic or mechanical means 
provides an efficient approach to engineer the band structure of 
MoS 2 . 
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